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Abstract 



A generic expression of double-parton scattering (DPS) cross sections in high-energy nucleus-nucleus collisions is 
derived as a function of the corresponding single-parton hard cross sections. We consider the case of prompt- //i/f 
production and find that the DPS contributions represent (15-20)% of the total JjiJ/ cross section in lead-lead (Pb- 
C ' Pb) collisions at the CERN Large Hadron Collider In Pb-Pb at 5.5 TeV, about 300 double-7/i^ events are expected 

per unit midrapidity in the dilepton decay modes, per inverse-nanobam. Double-y/i/^ DPS production can partially 
explain the relatively less suppressed yields of JjiJ/ observed in Pb-Pb collisions at the LHC compared to results at 
lower collider energies. 



, Introduction - The production of heavy-quark bound states of the charmonium and bottomonium (T) 

famiUes in high-energy proton-proton (p-p) and nucleus-nucleus (A- A) collisions is governed by both perturbative 
Q^l and non-perturbative aspects of Quantum Chromodynamics (QCD) and has been extensively studied at fixed-target 

D ' and collider energies |[T]. A pair of charm or bottom quarks (cc, bb) is first produced mostly in a hard gluon-gluon 

collision with cross sections computable via perturbative QCD (pQCD) calculations. The subsequent evolution 
of the QQ pair towards a color-singlet bound state is a non-perturbative process described in various theoretical 
approaches including color-singlet and color-octet mechanisms, non-relativistic QCD effective field theory, or color 
' evaporation models (see e.g. 121 for a review). 

in ■ 

, In the case of A-A collisions, quarkonium has been since long proposed as a sensitive probe of the thermo- 

dynamical properties of the hot QCD medium produced in the course of the collision [iSJ. Analysis of quarkonia 
correlators and potentials in finite-temperature lattice QCD l?) indicate that the different cc and bb bound-states dis- 
sociate at temperatures T for which the color (Debye) screening radius of the medium falls below their correspond- 
ing QQ binding radius. Experimental confirmation of such a quarkonia dissociation pattern would provide a direct 
means to determine the QCD transition temperature possibly reached in the system [jSj. Somehow surprisingly, 
J I if/ production in lead-lead (Pb-Pb) collisions at the LHC is observed to be less suppressed - compared to 
baseline p-p collisions at the same energy - than at the Relativistic Heavy-Ion Collider (RHIC) [T0| despite the fact 
, that the average medium temperature at LHC nucleon-nucleon center-of-mass (cm.) energies ( ^fs^ - 2.76 TeV) is 

^ ■ at least 30% higher than at RHIC ( y[s^ - 200 GeV) (TTl. At the LHC, the Jjifj suppression factor at low transverse 

- - -' momenta pj amounts to 7?aa ~ 0.5 whereas it was measured to be as low as /?aa ~ 0.2-0.3 at RHIC. Approaches 

combining Jjtp dissociation in a deconfined phase plus regeneration due to charm-quark recombination llT2l can 
reproduce the observed trends in the data although the model parameters {ctcc cross section, medium density, ...) 
need to be validated with other LHC observations. 

In this paper we discuss and quantify the role of double-parton-scattering (DPS) processes as an extra source of 
J 1 4/ production in Pb-Pb collisions at LHC energies. Due to the fast increase of the parton flux at small parton frac- 
tional momenta, x = Pparton / /'hadron lfT3l , the probability of having multiple hard parton interactions (MPI) occurring 
simultaneously at different impact parameters increases rapidly with collision energy and constitutes a significant 
source of particle production in p-p |1I4] and, in particular, nucleus-nucleus [15] collisions at semihard scales of a 
few GeV. The evidence for double parton scattering (DPS) processes producing two independently-identified hard 
particles in the same collision is currently based on p-p and p-p measurements with final-states containing (i) multi- 
jets at V* = 63 GeV lO, 630 GeV UJ], and L8 TeV El, (ii) 7-i-3-jets at V* = 1-8 TeV |[19| and L96 TeV 1201; 
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(a) (b) (c) 

Figure 1: Schematic DPS contributions in A-A collisions: (a) The two colliding partons belong to the same pair 
of nucleons, (b) partons from one nucleon in one nucleus collide with partons from two different nucleons in the 
other nucleus, and (c) the two colliding partons belong to two different nucleons from both nuclei. 



and (iii) W+2-jets at ^fs - 7 TeV ||2T]| . Such measurements show an excess of events in various differential dis- 
tributions with respect to the expectations from contributions from single-parton scatterings (SPS) alone. LHC p-p 
measurements of double-// 1/' production |'22j as well as of single-7/i/' production as a function of the event multi- 
plicity II23I have been also interpreted in the context of DPS ll24H27l and MPI models respectively. The 0(10%) 
"Cronin" enhancement observed in inclusive hadron production above pr ~2 GeV/c in p-Pb at the LHC ||28]| can be 
also naturally explained as due to extra DPS contributions (although, the fragmentation of the partons into final-state 
leading hadrons, carrying out a fraction of the original momentum, washes out the effect). 

We investigate DPS in A-A collisions following a previous similar study for proton-nucleus collisions ||291 . The 
larger transverse parton density in nuclei compared to protons results in enhanced A-A DPS contributions coming 
from interactions where the two partons belong or not to the same pair of nucleons of the colliding nuclei (Fig.[Tli. 
Consequently, in Pb-Pb at LHC energies we expect a non-negligible probability of two parton-parton interactions 
independently producing two J /if/ mesons in the same nuclear collision. 

Cross sections for double parton scattering in proton-proton and nuclear collisions - The DPS cross section 
in p-p collisions can be obtained from the convolution of parton distribution functions (PDF) and elementary cross 
section summed over all involved partons (see e.g. Il30l ) 

= (f )Z [ Oxi^-2;bi,b2;e?,e^) 

ij,kj 

x&%,,x\,Q\)&l{x2,x'^,Ql) (1) 
X T"{x[,x'2, bi - b, b2 - b; Q\, Q\) dxidx2dx[dx'2d^bid^b2d^b , 

where rp(xi, X2; bi, b2; Q^, are double-PDF which depend on the longitudinal parton momentum fractions x\ 
and X2, and on the transverse position bi and bj of the two partons undergoing the hard processes at the scales 
Q\ and Q2, 0-'^ and are the parton-level subprocess cross sections, and b is the impact parameter vector con- 
necting the centers of the colliding protons in the transverse plane. The combinatorial factor m/2 accounts for 
indistinguishable (m = 1) and distinguishable (m - 2) final-states. Decomposing the double PDF into longi- 
tudinal and transverse components, with the transverse components expressed in terms of the overlap function 
/(b) = J /(bi)/(bi - h)d^bi for a given parton transverse thickness function /(b), and making the assumption that 
the longitudinal components reduce to the "diagonal" product of two independent single-PDF, the cross section of 
double parton scattering can be finally expressed in the simple generic form IIT4l 



where cr^^^ is the inclusive single-hard scattering cross section, computable perturbatively to a given order in as, 

^ J D'p(xi;Ql)f(bi)d-'!;(xux[) x o'^ix'^; Ql)f(bi - b)dxidx[d^bid^b 

^ j D'pixi ; Ql) a-*(xi , x[) D\{x\ ; Q\)dx,dx\ , (3) 



- z 

and cr^i, is a normalization cross section representing the effective transverse overlap area of partonic interactions 
that produce the DPS process 



J d-bt- 



(b) 



-1 



(13 + 2) mb, (4) 



whose numerical value has been obtained empirically from fits to p-p and p-p data 1(161 - 1211 . 

To compute the DPS cross section in nucleus-nucleus collisions we proceed as done for proton-nucleus colli- 
sions in 1291 . The parton flux is enhanced by the number A of nucleons in each nucleus and - modulo (anti)shadowing 
effects in the nuclear PDF (see below) - the single-parton cross section is simply expected to be that of p-p colli- 
sions, or more exactly that of nucleon-nucleon collisions (N-N), scaled by the factor A~, i.e. 



SPS ^ SPS 

{AA^ab) " (NN- 



.ah) / TA(bi)TA(bi - b)d2b,d2b = trfP^^^,, J TAA(b)d2b = • a-f^l^.^,^ . (5) 



Here rA(b) is the nuclear thickness function as a function of the impact parameter vector b connecting the centers 
of the colliding nucleus in the transverse plane, and TAA(b) the standard nuclear overlap function normalised to 
II3TI . The DPS A-A cross section is thus the sum of three terms, corresponding to the diagrams of Fig. ((T); 

1. The first term corresponding to Fig. ^a), similarly to the SPS cross sections Eq. (|5]l, is just the DPS cross 
section in N-N collisions scaled by A~ 

^{AA^ab) - ^ ■ ^NN^ah) ■ 

2. The second term. Fig. [Ttb), accounts for interactions with partons from one nucleon in one nucleus with 
partons from two different nucleons in the other nucleu^ , 

DPS,2 _ ^ DPS X /-7\ 

'^(AA^ab) - -^^(NN^ab) ' ^.ir.w ' ^ 2.AA > (') 

with 

T2,AA = J TA(bi)TA(bi - b)TA(bi - b)d2bid2b = (A - 1) J d2rTi(r) - (A - 1) ■ Taa(O). (8) 

3. The third contribution from interactions of partons from two different nucleon in one nucleus with partons 
from two different nucleons in the other nucleus. Fig. [TJc), reads 

DPS, 3 _ J)PS /QN 

"(AA^ab) ~ "(NN^ab) ' "cs.pj, ' ^ 3,AA , K^) 

with 

T3,AA = (^^j jTA(bi)TA(b2)TA(bi-b)TA(b2-b)d2bid2b2d2b = J^^J ^A\Tl^{r). (10) 

where the integral of the nuclear overlap function squared does not depend much on the precise shape of the 
transverse parton density in the nucleus and it amounts to .94 ■ Taa(O) for the hard-sphere approximation 
and A^/2 ■ Taa(O) for a Gaussian profile. 



'Possible extra "non-diagonal" interference terms, computed for light nuclei in 1321 , would lead to an increase of the value of o"^^^^ but 
are effectively covered by the quoted uncertainty for (T^^^^, Eq. j4). 
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The factors (A - 1)/A and [(A - 1)/A]^ in the two last terms take into account the difference between the number 
of nucleon pairs and the number of different nucleon pairs. Adding (|6]l, (|7]i and (|9]l, the inclusive cross section of a 
DPS process with two hard parton subprocesses a and b in A-A collisions (with A large, so that A - 1 ^ A) can be 
written as 



DPS 
" (AA^ab 



— /I cr,. 



PS 



(NN^ab) 



A 2 DPS 

" " (NN^ab) 



, , 2(A-1) 

1 + ^2 ^..r,p 



/ d^rriCr) + (^j cr^,^^ J d2rTi^(r) 



A ' 



,,pTaa(0) + -(r,,„TAA(0) 



(11) 



In the simplest hard-sphere approximation for a nucleus with uniform nucleon density of radius Ra = r^A^^^ and 
ro = 1.25 fm, the nuclear overlap function at ^ = amounts to Taa(O) = 9AV(87rR2) = 31.5 mb^' for ^o^pb -^"Spb 
collisions. A direct evaluation of the integral using the measured Fermi-Dirac spatial density for the lead nucleus 
(Ra = 6.624 fm and surface thickness a = 0.546 fm) ^ yields Taa(O) = 30.4 mb^'. Using the latter Taa(O) 
value and cr = 13 mb, we compute the expression in parentheses in Eq. (flU - which quantifies the total DPS 
enhancement factor in A-A compared to N-N collisions, Eq. (|6]l - to be of the order of 200, dominated by the hard 
double nucleon scattering contributions. Fig. [Tic). The final DPS cross section "pocket formula" in nucleus-nucleus 
collisions as a function of the elementary nucleon-nucleon single-parton cross sections, can be obtained combining 
Eqs. dill and ( fTTT i: 



cr, 



_DPS 

{AA- 



.ab) - ( 2 ) 



SPS . SPS 

" (NN^a) " (NN^b) 



(12) 



with the effective A-A normalization cross section for Pb-Pb amounting to 



Creff.AA - 



A2[l 



, Taa(O) + 



■ 0"cir.pp TaA 



(0)] 



= (1.5 + 0.2) nb. 



(13) 



Numerically we see that whereas the single-parton cross sections in Pb-Pb collisions, Eq ((Sj, are enhanced by a 
factor of A^ ^ 4 ■ 10"* compared to that in p-p collisions, the corresponding double-parton cross sections are 
enhanced by a much higher factor of 0-,,^ /o-gj aa - 9 ■ 10^. 



Results - From Eqs. (fT2l i. with m = 1, and ( fT3] ) we can compute the expected double-parton cross sections for 
y/i^-pair production in Pb-Pb collisions using the single-parton J/tfr cross sections in nucleon-nucleon collisions, 
'^m^j/i/'X)' obtain the latter quantity at all relevant collision energies, we use the color evaporation model 
(CEM) predictions of |(34|, properly scaled to match the existing p-p and p-p data, and taking into account the 
shadowing of the nuclear PDF with respect to the free nucleons PDF ||35]| . The SPS cross sections for prompt J/t// 
have been measured down to zero pr in p-p at ^fs = 1 .96 TeV at rapidities |y| < 0.6 ll36l and in p-p at ^fs = 2.76 TeV 
(lyl < 0.9 El, 2 < lyl < 4.5 IM) and 7 TeV (|y| < 0.9 ||39l, 1.6 < \y\ < 2.4 gOl, 2 < lyl < 4.5 El]) after subtraction 
of the decay contributions from bottom meson^. The extrapolations to total J/t// cross sections at the LHC have 
been obtained by integrating a Gaussian distribution fitted to the data points measured at different y. The Tevatron 
midrapidity cross section has been extrapolated to full-y based on the prescription proposed in lB3l . The obtained 
values, with their propagated uncertainties, are listed in Table[T]and shown as data points in Fig.|2](top). 

The corresponding CEM predictions for crfpp^j/^x) scaled by a K-factor of twc0to match the experimental 
results, are also plotted as a solid curve. The values for cr^^l/^j/^x) obtained from the CEM predictions scaled 
to the p-p data and normalized to the ratio of the p-p over N-N cross-sections, where the latter have been determined 
using the EPS09 nuclear PDF 1451 . In the x ^ lO^^*, ^ mj^^ region relevant for J/tfr production at the LHC, the 
Pb gluon PDF is moderately depleted, by a factor of (1 -5g_pb) ~ 15%-25% with respect to the free nucleon density, 
resulting in a reduction of the gg J/ilr +X yields compared to the p-p data by a factor of ( 1 - S^p^) ~ 30%^0% 
(dashed-dotted line in Fig. |2l top). The EPS09 uncertainties, of the order of +10-15%, have been propagated in 
quadrature with those associated with the data extrapolation, into the N-N cross sections quoted in Table [1] 

^Decay contributions from higher-mass and i//' resonances increase the prompt- J/i^ yields by about 35% 1421 . 

^More recent CEM calculations |44| , with higher values of the c-quark mass and lower renormalization scales fig, reproduce well the LHC 
data. Note also that feed-down ifr' contributions are not included. 
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Table 1: Total cross sections at LHC energies for the production of prompt Jjtj/ in single-parton-scatterings (SPS) 
in p-p, N-N, Pb-Pb collisions, and of prompt y/i/r-pairs in double-parton-scatterings (DPS) in Pb-Pb. The p-p 
values are extrapolated from experimental data, the N-N values are based on the CEM (normalized to p-p,p-p 
data) including EPS09 nuclear PDFs, and the Pb-Pb results are derived from the N-N cross sections via the quoted 
equations. 



Process 


Cross section 






(TeV) 








1.96 


2.76 


5.5 


7.0 




measured (extrapolated) 


25. ± 9. 


28. ± 8. 




49. ± 9. 


o-SPS(N-N^ y/i/rX) [yub] 


CEM+EPS09 PDF, Eq. © 


15. ±5. 


19. ±6. 


26. + 8. 


29. ± 9. 


crSPS (Pb-Pb ^ //l/rX) [mb] 


Eq. © 


650 + 200 


850 ± 250 


1 100 ± 350 


1250 ± 350 


^DPS (pb.pb ^ JI^JI^X) [mb] 


Eqs. dH-dnil 


75 + 28 


120 ± 45 


225 ± 80 


280 ± 100 



^10" 




PbPb ^ J/i|/+X 


10^ 




PbPb ^ J/i|/J/i|/+X 


10^ 










10^ 










10^ 










10^ 




□ pp,pp J/\|/+X (data) 

pp J/v|/+X (CEM X k-f actor) 

- - NN^ J/\|/+X(CEMxk-factor, EPS09 PDF) 


10^ 








^ B 


10 

















^ 0.5| 1 

^ 0.4- -^--^ PbPb -^tJ/vJ/¥+X]/y 

io.3- 

03 

I 0.2 - ______ 

I °1 23456789 

is (TeV) 

Figure 2: Top: Prompt-7/ i/r production cross sections in single-parton interactions in p-p, N-N, and Pb-Pb collisions, 
and for double-parton J/tfrJ/ip in Pb-Pb, as a function of cm. energy. Bottom: Fraction of prompt-y/i/r events 
where a pair of J/t// is produced in Pb-Pb collisions, as a function of cm. energy (the band indicates the size of the 
associated uncertainty). 
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The uppermost curves in the top panel of Fig.|2]show the resuhing Pb-Pb cross sections for single- and double- 
Jjij/ production, whereas their ratio is shown in the bottom panel. Single prompt- y/i/r cross sections amount to 
about 1 b at the LHC (dashed curve), and 10% to 30% of such Pb-Pb collisions are actually accompaniecQ by the 
production of a second J I if/ from a double-parton-scattering interaction (dotted curve). It is interesting to note that 
the rise of the DPS/SPS ratio appears to slow down at the highest ^fs^ as the shadowing of the nuclear gluon PDFs 
(which enters squared in the numerator but only linearly in the denominator) increase^ thereby reducing the total 
double-7/i/f yields. The uncertainty of the DPS fraction (shown as a yellow band in Fig.|2l bottom) amounts to about 
+20% including, in quadrature, the uncertainties from the EPS09 PDFs and from the effective cr^,f value, Eq. ( fT3l l. 
The dominant uncertainties in the absolute //i/' cross sections, coming from the experimental data extrapolation, 
cancel out in the single-to-double y/i/r ratio. 

The computed Pb-Pb cross sections are baseline theoretical results - for "minimum bias" (MB) collisions with- 
out any selection in the reaction centrality - including nuclear PDF modifications but no final-state effects (e.g. 
from color deconfinement and/or charm recombination) which can modify the final measured yields. Experimen- 
tally, Pb-Pb collisions at 2.76 TeV show a reduction factor of ~2 of the MB //i/r yields with respect to p-p, i.e. 
^AA ~ o-aa/(^^ ■ o"pp) ~ 0.5, whereas the corresponding value at RHIC amounts to x 0.4. At the LHC energy, 
according to Fig.|2](bottom), we expect a ~15% increase of the prompt J/if/ yields due to DPS contributions (which 
are negligible at RHIC), and thus for the same medium-induced prompt- //t/r suppression at both energies, the /?aa 
can be higher at the LHC due to the extra initial-state J /if/ production. The effect of the DPS contributions on the 
i^AA ratios is much more significant for central A- A collisions since the dominant component scales as the square 
of the nuclear overlap function, Eq. (fTOl i. The production of double-//)/' is thus strongly biased towards the most 
central collisions with lowest impact parameter This fact can also explain the much less depleted J/iJ/ yield in 
central Pb-Pb at 2.76 TeV (R"^^^ ~ 0.5) than at RHIC {R"^^ ^ 0.2-0.3). A quantitative analysis of the centraUty 
dependence of the DPS J/iJ/ yields is object of a coming study PTl . 

A confirmation of our predictions would be possible by measuring the cross sections for the concurrent produc- 
tion of two J ftp mesons in the same Pb-Pb interaction via their visible dilepton decay channels. At LHC energies, 
the cross section per unit-rapidity for single- //(/r amounts to dcrj/^/dy ~ at the (low-pr) rapidities covered 

by ALICE (at y - 0) and CMS (at y - 2), the detector acceptance and reconstruction efficiencies reduce the mea- 
sured yield by factors of ~12-14 OIlOl, and the dilepton branching ratio amounts to 0.0594. Squaring all these 
quantities for the case of double-7/ tfr production results in a final reduction factor of order 3 ■ 1 0"^ for both rapidity 
ranges. Thus, at 5.5 TeV one would expect a visible DPS cross section of about da^J'^j^^/dy\y=o(2) ~ 75 nb per 
dilepton decay mode, i.e. about 300 double-y/i/r events per unit-rapidity in the four combinations of dielectron and 
dimuon channels in 1 nb"' of integrated luminosity, assuming no in-medium suppression. The same estimates for 
the 15 (150) /.<b ' of Pb-Pb data already collected at 2.76 TeV result in about 2 (20) double-7/i/r events in ALICE 
(CMS) at mid (forward) rapidities. The very large combinatorial background of dilepton pairs with invariant masses 
around m j/^ makes however this measurement very challenging on an event-by-event basis (except maybe at higher 
Pt values where the background is smaller). 

Conclusions - We have derived the general expression for double-parton-scattering (DPS) cross sections in 
heavy-ion collisions as a function of (i) the elementary single-parton cross sections in nucleon-nucleon collisions, 
and (ii) an eff'ective ct^j^a parameter describing the transverse density of partons in the system. The DPS cross sec- 
tions in Pb-Pb are found to be enhanced by a factor of 200 compared to the A^-scaling of single-parton scatterings. 
We have studied the case of double-//)/' production at LHC energies and found that the DPS contributions represent 
(15-20)% of the total pmmpt-J/tp cross section and could well explain the relatively less suppressed J/if/ yields 
measured at ^/s^ = 2.76 TeV compared to results at lower collider energies. At 5.5 TeV, about 300 double-y/i/r 
events are expected per unit rapidity in the dilepton decay channels (in the absence of final-state suppression) for 
an integrated luminosity of 1 nb"', providing a quantitative test of the predictions presented here. 

''dps cross sections have to be understood as exclusive values which are part of the total inclusive single-parton cross section 1461 . 
^We note that both MPI and gg fusion (leading to shadowing of nuclear PDF) are counterbalancing processes characteiistic of the multi- 
scattering dynamics dominating at small- j:. 
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